Discrete-symmetry vortices as angular Bloch modes 
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The most general form for symmetric modes of nonlinear discrete-symmetry systems with non- 
linearity depending on the modulus of the field is presented. Vortex solutions are demonstrated 
to behave as Bloch modes characterized by an angular Bloch momentum associated to a periodic 
variable, periodicity being fixed by the order of discrete point-symmetry of the system. The concept 
of angular Bloch momentum is thus introduced to generalize the usual definition of angular momen- 
tum to cases where 0(2)-symmetry no longer holds. The conservation of angular Bloch momentum 
during propagation is demonstrated. 



Rotational symmetry in the x — y plane, denned by the 
0(2) symmetry group, implies the conservation of the 
z-component of angular momentum in two-dimensional 
nonlinear systems. If the 0(2) continuous symmetry is 
substituted by a discrete rotational symmetry, the usual 
approach based on Noether's theorem fails since infinites- 
imal symmetry transformations are no longer allowed. 
Conservation of angular momentum is expected to break- 
down in its usual form and, consequently, nonlinear dy- 
namics should show special features as compared to the 
0(2)-symmetric case. Since angular momentum is ex- 
pected to realize differently in discrete-symmetry sys- 
tems, it is particularly interesting to compare the differ- 
ent behavior of 0(2)-symmetric and discrete-symmetry 
vortices. A general feature of vortex solutions is their 
characteristic phase dislocation, which is determined by 
an integer number that will be referred to as vorticity 
(also known as winding-number, "topological charge" or 
spin). In 0(2)-symmetric systems, rotationally invariant 
vortices (i.e., vortex solutions whose amplitude is 0(2)- 
symmetric) are eigenfunctions of the angular momentum 
operator with eigenvalue (angular momentum) given by 
vorticity. Oppositely, discrete-symmetry vortices have no 
well-defined angular momentum. Examples of both types 
of vortices can be found in optics as well as in Bosc- 
Einstein condensate (BEC) systems. Optical vortices 
with 0(2)-symmetric have been experimentally observed 
in homogeneous self-defocussing nonlinear Kerr media Q 
whereas discrete-symmetry optical vortices have been ob- 
served in optically-induced photonic lattices 0, Q ■ The 
influence of discrete symmetry in the features of optical 
vortices are strongly reflected in their angular proper- 
ties. The phase of a discrete-symmetry vortex presents, 
besides the typical linear angular dependence character- 
istic of 0(2)-symmetric vortices, an additional sinusoidal 
contribution completely fixed by the order of the discrete- 
symmetry of the system (the so-called group phase) 
Besides, these systems exhibit a vorticity cutoff equally 
determined by the order of the discrete-symmetry, in such 
a way that no vortices of arbitrary order are permitted 



FIG. 1: (a) Mapping of the angular variable 9 into the non- 
compact coordinate s. (b) Equivalent periodic potential using 
this mapping (for an operator of the type L(s) = V 2 — V(s)) 
corresponding to a system with point-symmetry of order n = 
6. 



I will propose a new approach to the search for vortex 
solutions of a general nonlinear equation that describes 
not only the conventional optical and BEC cases but also 
more general situations. I will prove a general expression 
for symmetric vortex solutions by showing that they are 
angular Bloch modes. On the second hand, I will extend 
the conservation of angular momentum to the discrete- 
symmetry case by introducing the concept of angular 
Bloch momentum and by demonstrating its conservation 
during propagation. 

Let us consider the 2D nonlinear eigenvalue equation 
for a system possessing discrete point-symmetry of the 
Cn type (discrete rotations of order n) : 
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The aim of this paper is two-fold. On the one hand, 



where Lq stands for the field-independent part of the 
full differential operator (linear part) and Lnl(|0|) for 
the nonlinear part. There is no restriction about the 
specific form of Lq and Lnl- Usually, Lq depends on 
the transverse position coordinates through gradient op- 
erators and explicit functions of position (a typical ex- 
ample is Lq = V 2 — Vo(x,y)). The only assumption 
made is that the nonlinear part Lnl depends on the field 
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through its modulus exclusively. There can be also a 
explicit dependence on coordinates through other func- 
tions (e.g., inhomogcncous nonlinear coefficients: Lnl — 
72(2:, y)\(j)\ 2 + 74 (cc, y)\4>\ 4 ). It is considered that the sys- 
tem is invariant under discrete rotations C n , so that all 
operators and coefficients defining Lq and Lnl are as- 
sumed to be invariant under discrete rotations of nth or- 
der. If one is interested in symmetric solutions, i.e., those 
functions verifying that their modulus is a C ra -invariant 
(\4>(r,9 + 2ir/n)\ — \<fi(r,9)\), then the entire operator 
L = Lq + Lnl becomes invariant under discrete rota- 
tions: L(r, 9) = L(r, 6 + e), e = 2ir jn. One can transform 
this angular dependence into a dependence on a length 
variable — defined within an interval of the real axis — by 
using the mapping 6 — ► s = 9D/2n. This mapping maps 
the unit circle Si (where 9 is defined: 9 6 [— it, tt]) onto 
the real axis interval [-D/2, D, /2) (see Fig0. In terms 
of the new length variable s, invariance under discrete ro- 
tations of the L-operator becomes a periodicity property: 
if L(r, s) = L(r, 2tts/D) then L(r, s + a) = L(r, s) where 
a is the period in the non-compact variable s. Accord- 
ing to the previous mapping, the period a is fixed by the 
discrete-rotation angle e and thus it depends on the or- 
der of the discrete symmetry: a = cD/2tt = D/n. In this 
way, Eq.Q for symmetric functions is transformed into 
L(r, s)(p(r, s) — —£ip(r,s) where <p(r,s) = <fi(r, 2ns/D). 
One has thus converted the original discrete-symmetry 
eigenvalue problem into a problem of finding the self- 
consistent modes of a periodic operator. The modulus 
of these modes are invariant under periodic translations 
(\<p(r, s + a) I = \p{r, s)\), a consequence of the invariance 
of the modulus of a symmetric solution under discrete- 
rotations. As a clarifying example, one can consider an 
standard operator of the type L(r,s) = V 2 — V(r,s), 
where V(r,s) = Vo(r, s) + V^i l (\<p(i', s)\). If <p>{r,s) is a 
solution of Lip = —Sip with periodic modulus then it has 
to be en eigenmode of the periodic operator generated by 
itself L(r, s) = V 2 — V(r, s). According to Bloch's theo- 
rem, since V is periodic V(r, s + a) = V(r, s) the solution 
of an eigenvalue equation of the previous type is given by 
Bloch modes in the periodic variable s. This argument 
equally applies to the general periodic operator L(r,s), 
so that the most general solution of Lip = —Sip is given 
by ID Bloch modes 0: ip vv (r,s) = e tps u pl/ (r, s), where 
p is the Bloch momentum or pseudo-momentum of the 
mode and u p is the so-called Bloch function. The index 
v is known as the band index and occurs because, at a 
given value of p, there are many different eigenmodes of 
L. Since v does not play a role in the current discus- 
sion, I will omit it from the notation (although one has 
to recall that it is always present). The Bloch function 
is a periodic function of s with the periodicity of the L 
operator : u p (r, s + a) = u p (r, s). 

One can obtain interesting properties of solutions of 
Eq.Q by re- interpreting well-known properties of ID 
Bloch modes. A Bloch mode is characterized by the value 



of its Bloch momentum p. However, unlike plane waves, 
p is not an eigenvalue of the momentum operator and its 
value is constrained to lie in a restricted interval, called 
the Brillouin zone, defined by the condition \p\ < n/a. 
Moreover, due to the definition of the length coordinate 
s in terms of the angular variable 9, the ip p function has 
to be additionally periodic in the interval length D (s S 
[—D/2,D/2]),<p p (r,s) — p p (r,s + D). As a consequence, 
a discretization condition for the Bloch momentum p is 
obtained: e lpD = 1 => p = p m = 2nm/D, m E Z. If 
one reverses the 9 — > s mapping by re-introducing the 
angular variable 9 in p> p (r,s) = e tps u p (r, s), one obtains 
the solution in its original form: 

ct> m (r,9) =e ime u m (r,0), m e Z, \m\ < n/2. (2) 

I shall call these solutions angular Bloch modes. They 
constitute the general symmetric solutions of the non- 
linear eigenvalue equation Because of the angular 
nature of the 9 variable and its relation to the Bloch 
momentum p rn — 2irm/D, the index m will be referred 
to as the angular Bloch momentum (or pseudo angular- 
momentum) of the angular Bloch mode. The restriction 
on the permitted values of the angular Bloch momen- 
tum m is a consequence of the Brillouin zone limita- 
tion. The existence of the Brillouin zone for the Bloch 
momentum p establishes a condition for its permitted 
values: \p\ < n/a. Since p is discretized according to 
p rn — 2irm/D, m G Z, the Brillouin zone limitation effec- 
tively imposes a strict constraint into the angular Bloch 
momentum \m\ < D/2a. Inasmuch as the period a is 
determined by the order of symmetry n (a = D/n), one 
finds that the modulus of the angular Bloch momentum 
presents the upper bound \m\ < n/2 occurring in Eq.pjl. 

The expression J2J is the same one that it is ob- 
tained for vortex solutions in a nonlinear system with 
point-symmetry C n using group theory arguments 0, |]| : 

<t>i(r, 9) — e^ S (pQ {r,9), T S Z being the index of the group 

representation where the solution belongs to and (j)^ a 

discrete- rotation invariant function, </>q (r, 6 + 2n/n) = 

4>o(r, 9). The invariant function <j$ predicted by group 
theory is nothing but the Bloch function u m of the an- 
gular Bloch mode J2J|. Vorticity equals the index of the 
representation, so that 4>i represents a vortex with vor- 
ticity v = I (the exception is the I = n/2 — even n — 

solution) 0- The comparison of <pf(r,d) — e^ e <f>Q (r, 6) 
and Eq.(j2J) yields to interesting equivalences. Vortex soli- 
tons of vorticity v — I appear now as angular Bloch 
modes carrying angular Bloch momentum m = v = I 
(rn = n/2 — even n — excluded). The restriction on 
angular Bloch momentum \m\ < n/2, along with the 
previously established relation between vorticity and an- 
gular Bloch momentum, fixes a cutoff for the vorticity: 
\v\ < n/2 (v — n/2 — even n — excluded). This condition 
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FIG. 2: Schematic representation of Brillouin zones for two 
different order of symmetries: (a) even n (n — 6) ; (b) odd n 
(n = 5). Permitted angular Bloch modes are shown. 



can be also rephrased as: 

H < n/2 (even n) and \v\ < (n - l)/2 (odd n). (3) 

The explicit distinction between even and odd orders is 
made to remark the different behavior existent in the 
corresponding Brillouin zones in both cases. For odd n, 
the discretized Bloch momentum p m can never achieved 
the zone limit \p m \ < ir/a, whereas this is certainly pos- 
sible for even n. This can be clearly visualized in Fig. 
121 where two examples of Brillouin zones for even and 
odd n are shown. It is remarkable that the vorticity cut- 
off expressed in Eq. J3J has been obtained in a completely 
different mathematical framework using group-theory ar- 
guments 

At this point, it is interesting to compare the discrete- 
symmetry situation with that corresponding to contin- 
uous rotational symmetry. In the 0(2) case, Noether's 
theorem ensures the conservation of the z- component of 
the angular momentum associated to the cj> field during 
propagation: 
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(4) 



This quantity is also the expectation value of the angu- 
lar momentum operator L z = —id/ 89 for the <fi field, i.e., 
jz = {<fi\Lz\4>)- Since 0(2)-symmetric vortex solutions are 
of the type <j>i(r, 9) = er fj(r), vortex solutions are eigen- 
functions of the angular momentum operator with eigen- 
value I. Therefore, j z equals / for (normalized to one) 
vortex solutions and both represent the same constant 
of motion. On the other hand, the vorticity v of 0(2)- 
symmetric vortices is directly given by I: v = j z = I. 
Thus, angular-momentum and vorticity are equivalent 
quantities for 0(2)-symmetric vortices. 

In a discrete-symmetry system the z-component of 
angular momentum j z (En.ljljl) is no longer conserved. 
Noether's theorem does not apply because the symmetry 
under consideration ceases from being continuous. It is 



then necessary to find a different conserved quantity that 
explicitly manifests the discrete-symmetry invariance of 
the system. It will be proved next that this quantity is 
nothing but the angular Bloch momentum. In order to 
carry out the previous demonstration, it is simpler to re- 
sort to the representation of Bloch modes in terms of the 
non-compact variable s. It will be assumed that evolution 
is provided by a first-order operator in the evolution vari- 
able z, i.e., Lcp(r,s,z) = —idtp(r,s,z)/dz. This type of 
evolution operator includes Schroedinger-like equations 
(like those appearing in nonlinear paraxial evolution in 
optical systems or in the evolution of Bose-Einstein con- 
densates) but also more general evolution equations, as 
non-paraxial forward equations in axially-invariant opti- 
cal systems @- My aim now is to analyze the evolution 
of a field amplitude whose initial condition is given by a 
Bloch mode with well-defined pscudo- momentum p, i.e., 
(p(r, s,0) = ip p (r, s) = e lps u p (r, s). It is important to 
stress that this mode does not need to be a stationary 
solution of the vortex type satisfying the nonlinear eigen- 
value Eq.JU, so that its amplitude can evolve with z. It is 
a known fact that Bloch momentum is conserved during 
propagation in a periodic linear system. However, this 
question is not obvious in the nonlinear case inasmuch 
as the evolution operator explicitly depends on the field 
itself. 

Let us consider the total evolution from to z as a 
succession of infinitesimal evolution steps of length e in 
the limit e — ► 0. In a first-order evolution equation, the 
infinitesimal evolution of a field amplitude from the axial 
slice Zj to the following slice Zj + \ = zj + e is given by: 



(5) 



Now I will prove that if tpo (ipo = tp(r, s, 0)) is a Bloch 
mode with well-defined Bloch momentum p, then ipj is 
also a Bloch mode with the same pseudomentum p for all 
values of j. The proof is carried out by induction. Let us 
start by calculating the pseudo-momentum of the j = 1 
amplitude tpi. The initial amplitude tpo is a Bloch mode 
with pseudo-momentum p (ipo — ip p (r,s) = e lps u p (r, s)) 
and thus it is an eigenfunction of the translation oper- 
ator T a : T a ip (r,s) = ip (r, s + a) = e ipa ip Q {r, s). The 
amplitude \ipo\ is a periodic function of s with period- 
icity a since it is the modulus of the periodic Bloch 
function u p (r,s): \<po{r, s + a) = \(po(r,s)\, so that the 
nonlineal part of the L operator is translational invari- 
ant, L^ h (\ip {r, s + a)\) = L NL (\ip {r, s)|). Therefore 
the full operator L is invariant under finite translations 
since both its linear and nonlinear part are invariant, 
|X(|<^o|); T a ] — 0. This implies in turn that the infinites- 
imal evolution operator in Eq.|0 commutes with T a : 
[ e i£(M)e iTo ] [1 + iL(\<po\)e,T a ] = 0. Now I ap- 

ply the T a operator onto ipi and take into account Eq.© 
for j = to find T a L Pl = e^^T^o = e ipa ipx, where 
I have used the fact that the evolution and translation 
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operator commute between them and that ipo is an eigen- 
function of T a . Consequently, (pi is a Bloch mode with 
pseudo-momentum p. Now, following the induction pro- 
cedure, I will assume the Bloch-mode property to be true 
for ipj and I will prove it for ifj+i- The demonstra- 
tion is analogous to the j — case. I assume ipj is a 
Bloch-mode with pseudo-momentum p. Then its modu- 
lus is translational invariant, so LQipj |) and e lL ^^ e are: 
[L(\ipj\),T a ] = [e iL ^^ E ,T a ] = 0. If one acts with the 
translation operator T a onto tpj+i, commutes the transla- 
tion and evolution operators and takes into account that 
T a (fj — e lpa tpj (Bloch-mode condition for ipj), one read- 
ily finds that T a tfj +1 = e ipa (pj+i. This shows that (fj 
is a Bloch-mode with pseudo-momentum p for all values 
of j if the initial amplitude tpo is. Expressed in differ- 
ent words: axial nonlinear evolution preserves the Bloch 
momentum p. 

The continuous limit s — > of the previous statement 
permits to give an expression for the evolving field am- 
plitude cp(r, s,z): if ip(r, s,0) = e lps u v {r 1 s) then it is 
also true that ip(r, s,z) — e tps u p (r, s, z), u p being a in- 
dependent Bloch function. If we reintroduce the angular 
variable 9 instead of the non-compact one s, one obtains 
that the general expression for the evolution of an ini- 
tial angular Bloch mode will well-defined angular Bloch 
momentum m in the system under consideration is: 

<l> m {r,e,z) = e im9 u m {r,6,z), meZ, \m\ < n/2. (6) 

This equation establishes that the angular Bloch momen- 
tum m is preserved by nonlinear evolution. 

The interplay between ordinary angular momentum 
and angular Bloch momentum can be visualized by cal- 
culating the expectation value of the angular momentum 
operator for the Bloch mode amplitude ©■ One easily 
finds that j z (z) = m + {u m \L z \u m ) (z) or 

m = iz{z) - (u m \L z \u m ) (z). (7) 

This property shows that angular Bloch momentum 
is conserved despite angular momentum is not. The 
expectation value appearing in Eq.Q is j™(z) = 
Jr2 u^ n (—id/d9)u m corresponding to the angular momen- 
tum associated to the ^-dependent Bloch functions u m . 
The conservation of the angular Bloch momentum is then 
the result of a subtle balance between two non-conserved 
quantities: the conventional angular momentum j z (z) 
and the angular momentum j™(z) related to Bloch func- 
tions. The latter can be attributed to the presence of the 
discrete-symmetry system that acts as an angular peri- 
odic crystal. The dependence of the u m Bloch functions 
on 9 is a consequence of the angular periodicity of the sys- 
tem. In order to see it, it is clarifying to consider an 0(2)- 
invariant medium as the n — > oo limit of a system with 
discrete symmetry of nth order. In such a case, the u m 
Bloch functions become independent of the 9 angle and, 



therefore, their associated angular momentum vanishes 
fn^) ~ * 0. As a consequence, angular momentum 
becomes angular Bloch momentum j z (z) m an d, 

therefore, a constant of motion. It is then reasonable to 
attribute the angular momentum contribution j" l (z) to 
the discrete-symmetry system acting as an angular peri- 
odic crystal. One can interpret Eq.lQ) as the statement 
that during propagation conventional angular momen- 
tum is transfered to the discrete-symmetry system (the 
angular crystal), and vice- versa, in such a way angular 
Bloch momentum is conserved: dj z (z)/dz = —dj™(z)/dz. 

All the previous arguments apply to every type 
of evolving field provided it satisfies the initial con- 
dition of having well-defined angular Bloch momen- 
tum. Conservation of angular Bloch momentum ap- 
pears then as a property which occurs independently 
whether the solution is stationary or not. Nonlinear sta- 
tionary solutions of Eq.([Q) are only particular cases to 
which the general conservation law can be applied to. 
For discrete-symmetry vortices the angular-momentum- 
vorticity equivalence v — j z of 0(2)-invariant solutions is 
lost. Since discrete-symmetry vortices are angular Bloch 
modes of vorticity v = m, Eq.lJjJ establishes a new rela- 
tion between vorticity and the angular momentum car- 
ried by the vortex field (as defined in Eq.QJ): v = j z — j™. 
In this way, the angular Bloch momentum concept per- 
mits to unveil the role played by angular momentum and 
vorticity in systems whose symmetry is no longer 0(2) 
but a point-symmetry of lesser order. 
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